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ABSTRACT 

In this paper we describe a set of models to predict the colours of galaxies over a 
wide range of redshifts. We present example output from the simulations, and discuss 
their application to the selection of galaxies at high redshifts, particularly through 
identification of the Lyman break. Additionally we consider the optimal choices of 
filters for selection at a range of redshifts. 

An interface to a subset of the simulations has been made available 
on the World Wide Web for the benefit of the community at the location 
http : / / www-astro . physics . ox . ac . uk/~re j s/research/ galcols . html. 

Key words: galaxies: formation - galaxies: distances and redshifts - galaxies: funda- 
mental parameters ~ galaxies: photometry 



1 INTRODUCTION 

In recent years, two-colour selection techniques have proved 
themselves to be an invaluable method for identifying galax- 
ies at high redshifts, paricularly by means of chosing a filter 
set spanning the wavelengths between Lya and the Lyman 
limit. Most notable has been the work of Steidel and collab- 
orators (e.g. Steidel et al. 1996; Steidel et al. 1999). 

As part of a project to identify radio-quiet companions 
to known radio galaxies in the redshift range 3 ^ z ^ 4.5, we 
required theoretical models to determine the optimal selec- 
tion criteria given our choice of filters. This paper discusses 
modelling of the expected colours for galaxies and the poten- 
tial for efficiently distinguishing them from foreground stars 
and galaxies through broad-band photometric techniques. 
The results of our observations are described in a compan- 
ion paper (Stevens, Lacy & Rawlings 2000). 

In our searches for Lyman-break galaxies (hereafter 
LBGs), we initially worked to the assumptions that their 
spectra will be approximately flat above the rest-frame Lya 
wavelength, somewhat diminished between Lya and the Ly- 
man limit, and have little or no flux below 912A. Such as- 
sumptions were used for our pilot study of a region of ap- 
proximately one square arcminute around the z = 3.80 radio 
galaxy 4C 41.17, described in Lacy & Rawlings (1996). 

Unfortunately such naiVe selection methods lend them- 
selves to significant contamination from other red objects 
such as galaxies at moderate redshift (for instance where 



the 4000A break is redshifted to the position of the Lyman 
break at the targetted redshift) or low-mass main-sequence 
stars. With the benefit of the models described in this paper 
and of Hubble Space Telescope images of the region stud- 
ied in our pilot study, which have subsequently been made 
available through the HST Archive, we conclude that the 
majority of the high-redshift candidates previously identi- 
fied are in fact foreground contaminants. 

For more efficient selection, a much more detailed anal- 
ysis will be required and the following points must be con- 
sidered within the models: 

(i) Intrinsic spectra of starbursts 

(ii) Attenuation due to intergalactic Hi absorption sys- 
tems 

(iii) Effects of dust extinction 

(iv) Broad-band colours and potential for confusion with 
other objects through photometric selection techniques 

(v) Apparent magnitudes of high-redshift star-forming 
galaxies for particular cosmologies, likely star-formation 
rates and consequences for their observability 

Throughout this paper, we follow the convention estab- 
fished by Steidel et al. (Steidel & Hamilton 1992) of using 
the AB magnitude system of (Oke & Gunn 1983), conver- 
sions for which are given in Table |l| This has the great ad- 
vantage of being a magnitude system in which the colours 
of an unreddened star-forming galaxy between bandpasses 
sampling the rest-frame ultraviolet continuum will be close 
to zero. 



* Present address: Oxford University Computing Services, 
Banbury Road, Oxford, OX2 6NN 
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Filter 


^eff (A) 


^AB (Vega) 


U 


3664 


0.780 


B 


4394 


-0.105 


G 


4869 


-0.098 


V 


5449 


0.000 


R 


6396 


0.186 


I 


8032 


0.444 



Table 1. Conversions between AB and Vega-based magnitudes 
of the principal filters used in our models. 

1.1 Previous work 

Guhathakurta and co-workers (Guhathakurta, Tyson & Ma- 
jewski 1990) searched for LBGs in relatively deep UBjRI 
images, but were unable to provide followup spectroscopic 
evidence. They estimated the strength of the break at the 
Lyman limit from consideration of the UV spectra of O and 
B stars, with synthetic spectra showing breaks of a factor in 
the range 2 — 5 for a wide variety of input initial mass func- 
tions (IMFs), even without considering interstellar or inter- 
galactic absorption by neutral hydrogen. Candidate LBGs 
are selected from their U — Bj and Bj — I colours. Their re- 
sulting selection region is somewhat conservative compared 
to later work, as it does not take into account the possibility 
of a substantial drop in flux below Lya. 

Most of the work of Steidel and collaborators makes use 
of a non-standard filter system, originally chosen for its suit- 
ability to a field containing a QSO and damped Lyman a 
system at known redshifts (Steidel & Hamilton 1992). Ini- 
tially, selection is done taking into account intergalactic ab- 
sorption measurements from QSO spectra (Sargent, Steidel 
& Boksenberg 1989), while later selection criteria are revised 
to take into account the detailed analysis of intergalactic at- 
tenuation by Madau (1995), which we use ourselves in § ^ 
Subsequent spectroscopic observations show that the selec- 
tion process is highly effective with less conservative selec- 
tion criteria. The lower bound in Un — G has been relaxed 
from 1.5 to 1.0, with a resulting increase in surface den- 
sity (to the magnitude limit of TZ < 25.0) from around 0.5 
candidate objects per arcmin^ to around 1.25 per arcmin^ 
(Steidel et al. 1997; Adelberger, Steidel et al. 1998). 

Madau et al. (1996) use a similar approach to that de- 
scribed here in order to determine appropriate regions in 
two-colour space for the selection of U300 and i345o dropouts 
in the northern HDF. We develop this approach to a wider 
range of filters, including those typical of those found on 
ground-based optical CCD and NIR cameras, and to span 
as wide a range of potential galaxy scenarios as possible. 



2 SIMULATED BROAD-BAND COLOURS OF 
GALAXIES 

In order to devise efficient selection criteria for high-redshift 
galaxies, it is necessary to model the expected colours of a 
wide variety of different types of galaxies for all redshifts at 
which one might possibly expect to be able to see objects. 

We base our models on the pegase code of Fioc & 
Rocca-Volmerange (1997) (hereafter FRV); our thanks go 



to them for making the code publicly available. [T| The FRV 
code produces synthetic galaxy spectra to a resolution of 
lOA from 220-8780A in the rest-frame|j and to 200A in the 
NIR. Like most current techniques for modelling the inte- 
grated spectrum of a population of stars, the FRV code uses 
the isochrone technique (Bruzual & Chariot 1991). A library 
of higher-resolution spectra is available, but this only covers 
the rest frame optical and is thus of little use when the ob- 
jects are redshifted; besides which the higher resolution data 
is not overly crucial when dealing with broad bandpasses. 

A substantial library of synthetic galaxy spectra using 
a wide variety of different input parameters to the FRV code 
was built up. 

2.1 Variations of input parameters 

2.1.1 Evolutionary tracks 

Two sets of stellar evolutionary tracks are made available 
with the FRV code, namely "Padova" (Bressan et al. 1993) 
and "Geneva" (Schaller et al. 1992; Charbonnel et al. 1996). 
Post main sequence evolution, in particular that along the 
asymptotic giant branch for the most massive stars (Jimenez 
et al. 1999), proves the most difficult phase to trace in such 
tracks. Primarily this affects the NIR luminosity, with only 
a slight effect in the UV. 

Another limitation is that the evolutionary tracks used 
are for solar metallicity, and will therefore not be an accurate 
representation of those in metal-poor regions such as may 
be found in the early universe. The version 2.0 FRV code at- 
tempts to redress this problem (Fioc & Rocca-Volmerange 
2000), and will be incorporated into our models at a future 
date. Early indications are that metallicity significantly af- 
fects the slope of the UV continuum, with bluer spectra at 
low metallicities. Evidently this will have a considerable ef- 
fect upon the derived star formation rates. Further work 
with the new models is required to investigate this more 
thoroughly. 

2.1.2 Initial Mass Function 

The FRV code by default incorporates the IMFs of Salpeter 
(1955), Scalo (1986), Miller & Scalo (1979) and of Kroupa 
(Kroupa, Tout & Gilmore 1993). The Salpeter is the simplest 
and still in widespread usage, but somewhat in error at the 
extremes. Later IMF models are an improvement in this 
respect, being somewhat more realistic and with less abrupt 
cut-offs at both ends, but are still limited by the the ability 
to observe the full range of stellar masses only in the solar 
neighbourhood. 

The Kroupa IMF results from work on stars in the 
galactic disk which resolves the apparent differences between 
observations made in the solar neighbourhood and those de- 
rived from surveys to much greater distances. It is probably 

t A more recent version of the code has recently been announced 
(Fioc & Rocca-Volmerange 2000), but at present our models are 
built around the version one code. 

i This resolution was higher than that of the Bruzual-Charlot 
models available at the time (Bruzual & Chariot 1993), an impor- 
tant factor when one wishes to shift the spectra to high redshifts. 
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Figure 1. 



Representations of tlie four initial mass functions used in the sim- 
ulations. For clarity, all have been multiplied up by a factor m^-^^; 
this gives a constant slope for the Salpeter IMF. 

to date the most realistic representation of the IMF in our 
own galaxy and it is used as the default in our models. 

Applying such empirically-derived IMFs to star- 
formation in the early universe is a big assumption, given 
our limited understanding of the physical processes involved. 
For example, low metallicity, different cloud dynamics and 
the possible influence of AGN activity upon star-formation 
may all affect the high-redshift IMF. 

The high-mass slope of the IMF is of particular inter- 
est to us since it determines the numbers of O and B stars 
which dominate the UV luminosity in star-forming galaxies. 
Any estimate of star-formation rate based on the UV lumi- 
nosity is thus strongly dependent on the choice of IMF. For 
simplicity and ease of comparison, star formation rates are 
generally derived for a Salpeter IMF, but one must bear in 
mind that with, for example, a Kroupa IMF as opposed to 
a Salpeter distribution, the overall star-formation rate may 
be a factor of '--^ 2 higher if derived from the UV luminos- 
ity. The effects of extinction due to dust will boost the SFR 
higher still. 

2.1.3 Nebular emission 

FRV use a default value for the fractional Lyman contin- 
uum absorption / = 0.7, in agreement with LMC Hll values 
(DeGioia-Eastwood 1992) , but for completeness values from 
to 1 were considered. This would have considerable effect 
on the strength of the Lyman break save for the fact that in- 
tergalactic attenuation becomes the dominant effect at red- 
shifts sufficiently high for us to observe below the Lyman 
limit (§ P). However it does govern the strength of emission 
lines, which, particularly in the case of Lya, can in principle 
have considerable effect on the observed colours, though it is 
often strongly suppressed. Continuum emission is increased 
in the rest-frame optical, but this will only have an effect on 
the infrared colours of galaxies at high redshift. 



2.2 Star-formation model 

A variety of star-formation models are included in the 
FRV code, ranging from instantaneous bursts to constant 
SFR scenarios. Primarily the choice affects how quickly the 
galaxy spectrum evolves from that of young starbursting 
regions to that of aged ellipticals. We include in the final 
models galaxy spectra covering a very wide range of star- 
formation models and ages, subject to a conservative cut- 
off to eliminate galaxies at ages unphysical for any reason- 
able cosmology given their redshift. Thus the models include 
more or less any conceivable galaxy type. 

2.3 Other parameters 

The models also incorporate variations of initial star for- 
mation rate, governing the timescale over which stars are 
formed, and the fraction of stellar ejecta available for fur- 
ther star formation. This latter only affects galaxies suffi- 
ciently aged for a substantial proportion of the stars origi- 
nally formed to go through the late stages of stellar evolu- 
tion; significant star formation resulting from stellar ejecta 
results in the overall galaxy spectrum being a combination of 
the UV spectrum of a young star-forming galaxy together 
with the optical and NIR spectrum expected of an aged 
galaxy. 

Dust extinction is also considered, but rather than than 
use the FRV extinction routine, we chose to use our own as 
described in § ^. 

The spectra were then convolved with selected broad- 
band filter profiles (those of the WHT optical and NIR filter 
sets, and the HST WFPC2 fifiers as used in the HDF obser- 
vations), scaled to the luminosity of a galaxy of stellar mass 
1O^"M0, and the AB magnitudes output. 



3 INTERGALACTIC ATTENUATION BY THE 
LYMAN-a FOREST AND BY LYMAN LIMIT 
SYSTEMS 

To approximate intergalactic attenuation, the method of 
Madau (1995) was followed. Equations 12 and 13 of his pa- 
per were used to correct for the Lyman- a forest. To correct 
for absorption of the Lyman continuum we made use of the 
approximate numerical integration of his Equation 16, which 
he states is good to 5%. For our purposes this is perfectly 
adequate and computationally inexpensive. 

Given the random distribution of the absorbers, signif- 
icant deviations from the mean attenuation are to be ex- 
pected at any particular wavelength. In the Lyman-a forest 
region, this does not present any problem, since the use of 
broad bandpasses means that one will be averaging over a 
large number of absorbers. 

Madau's figure 2 shows the measured flux decrements 
for quasar spectra within bandpasses of ~100A in the rest- 
frame. Few decrements differ from the estimated mean by 
more than about 0.1 mag. In our case we can expect an 
even greater averaging effect to act in our favour owing to 
the wider bandpass of our filters. 

Below the Lyman limit, a substantial contribution to 
the attenuation will come from the much rarer Lyman limit 
systems (LLS), for which the mean number along any line 
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Model spectra in fi, for a z = 3.5 galaxy undergoing star forma- 
tion at a constant rate of IOMq yr~^, and for an evolved galaxy 
at z = 0.5, typical of lower-redshift galaxies likely to have similar 
colours to the targetted high-rcdshift objects owing to the pres- 
ence of the 4000A break. Our selection criteria must therefore be 
able to discriminate effectively between the two, using the much 
bluer continuum colour of the high-redshift object relative to the 
low-redshift galaxy. The Lyman-a and Lyman-limit breaks are 
clearly visible in the former; similar objects may possibly pass 
the selection criteria for both U and B dropouts. 
Also shown are the response curves for various filters. The optical 
filters are those from the Isaac Newton Group telescopes on La 
Palma and include those used in our observations (a similar set of 
filters plus the Z is being used for the wide field survey of Dalton 
et al. currently in progress at the INT); all filter profiles have 
been convolved with the response curve for the TEK CCD used 
in our observations. The JHK filter profiles are taken from data 
available on the UKIRT WWW pages. 

No scan was available for the G filter (loaned by Richard McMa- 
hon); for these models a bandpass of width lOOOA centred at 
4900A has been assumed, as per the design specifications. 



of sight will be of order unity. There is, however, a signif- 
icant probability that there will not be a single LLS along 
a particular line of sight, about 10% at the redshifts of in- 
terest. While in the case of the Lya forest, the distribution 
of the number densities of absorbers about the mean can be 
approximated by a Gaussian at the redshift of interest, such 
an approximation cannot be made for the LLS. 

The chances of an object missing our selection cuts ow- 
ing to the absence of a LLS along the line of sight are greater 
for U dropouts than for B dropouts for two reasons. Firstly, 
the relative contribution of the LLS to the attenuation at 
2 ~ 3 is greater than at z ~ 4 (see Figure ^). Secondly, 
the expected number of LLS along a line of sight increases 
significantly between z ~ 3 and z ~ 4. 

M0ller & Jakobsen (1990) give values ranging from 0.25 
at Ze = 2.5 to about 0.05 at Ze — 4.5 for the probabil- 
ity of encountering no LLS leading to an absorption depth 
T > Tcrit ~ 1.5 at the Hen line Ae — 304A. As the photoion- 
isation cross-section of H decreases approximately as u'^ 
below the Lyman limit, this optical depth corresponds to a 
very optically-thick Lyman limit system, but these probabil- 
ities are a useful indicator of how the likelihood of encoun- 
tering a Lyman limit system changes with redshift. 
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Figure 3. 

Variation of magnitude decrements with redshift for the six broad- 
band filters used in our own observations. Solid lines show the 
mean decrement due to all absorbers, while dashed lines do not 
include the contribution of Lyman limit systems. Note in partic- 
ular that the absence of LLS absorbers along a particular line of 
sight is most likely to affect the chances of detecting U dropouts 
at redshifts 3 S: z S: 3.5. Even so, a substantial break should be 
observed even without a LLS absorber on the line of sight. 



4 DUST EXTINCTION 

While in the optical a measure of the dust extinction is ob- 
tainable from the emission line Balmer optical depth Tg, 
dust extinction in the UV is typically parametrized in terms 
of the continuum slope. The UV continuum of a star-forming 
galaxy may be approximated by a power law F{X) oc A", 
with a "colour excess" defined by i5(/3) — P — fio. Spectral 
synthesis models give a value —2.5 H /3o ~ —2.0, in broad 
agreement with measured values. For a young stellar pop- 
ulation (age ;S 2 X lO'^ yr), the differences, other than of 
normalisation, between the spectra of the extreme scenarios 
of instantaneous burst and constant star formation rate are 
small and within measurement uncertainties. 

We calculate dust extinction according to the method of 
Calzetti, Kinney & Storchi-Bergmann (1994, hereafter C94), 
taking the extinction parameter A{\) to be empirically fitted 
by a third-order polynomial Qe{x) (where a; = 1/A such that 
A{\) = t'bQ.{x). 

Qe is purely the selective extinction as opposed to the 
total extinction, since for any spectral slope (3 one does not 
in general know at what wavelength there is no extinction. 
To obtain the total extinction, we use the method described 
by Meurer et al. (1995), using consideration of energy con- 
servation (relating the far IR excess to /3). The ratio of total 
to selective extinction is parametrized as X'^ = Ax/E{P). 
For the HST F220W filter (central wavelength 2280A, thus 
Qe(F220W) = 1.582) and the C94 starbursts, this value is 
constant at X'p22ow ~ 



Using the starburst values of C94 to give 
0.494i5(/3) and choosing /3o = —2.5, we obtain 



/'(A) 



■j^Q-0.4xl.6(/3+2.5) jj-^-) 



exp [0.494(/3 + 2.5) (Qe(A) - 1.582)] 



(1) 



The dust extinction corrections are applied to our syn- 
thetic spectra prior to correcting for intergalactic extinction, 
using values for (3 of —2.5 (zero extinction), —1.5 and —0.5 
(equivalent to E{B - V) of 0.00, 0.23 and 0.46). 
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Figure 4. 

Variation of Z magnitude with redshift for an unreddened star- 
forming galaxy with a continuous star formation rate of 10 h^g 
Mq yr~^. The top plot shows the effect of changes in cosmology 
for a Kroupa IMF. The bottom plot shows the effect of varying the 
initial mass function for Q = 1, A = and Padova evolutionary 
tracks except where indicated. Sensitivity to the choice of evolu- 
tionary track is at the 5% level (though this is likely to increase if 
metallicity effects are taken into account), but the choice of IMF 
has a very significant effect on the luminosity-SFR calibration. 
Z is used as Lya does not reach the band until z fa 5.8, and 
hence it represents the unattenuated UV continuum flux. While 
the unreddened UV continuum is frequently regarded as a flat 
spectrum, this is merely an approximation, and K correction ef- 
fects lead to the variation of magnitude with redshift differing 
slightly from the distance modulus. 

5 LUMINOSITY-SFR CALIBRATION 

For the calibration between observed magnitude and star- 
formation rate, we choose for ease of comparison to use the 
same model as Steidel et al. (1996), namely a Salpeter IMF 
with upper mass cutoff at SOAf© — this low cutoff has a 
significant effect on the magnitude (Figure^. 

For an — 1, A = cosmology with Ho = 50 
km s~^ Mpc~^, we find that a constant star formation rate 
of IO.OM0 yr^^ with solar metallicity at z = 3.25 corre- 
sponds to a magnitude in R of 24.88 at age lOOMyr, tending 
asymptotically to J? = 24.79 as the galaxy ages^. This is 
slightly fainter than the value quoted by Steidel et al., but 
probably insignificant in view of the use of different spectral 
synthesis codes and stellar libraries. 

Figure ^ shows the effect of changes in IMF, cosmology 
and evolutionary tracks. It is evident that accurate deter- 
mination of the shape of the Madau curve is dependent on 
the cosmology, initial mass function and more accurate rep- 
resentations of stellar evolution at low metallicities. 



§ By an age of 60Myr, the luminosity is within 10% of this figure, 
and within 1% by an age of 400Myr. 




Figure 5. 

Variation of U — B and B — R colours with redshift z. Galax- 
ies are indicated by small dots coloured by redshift; young (age 
<100Myr) galaxies in the redshift range 2.75 < z < 3.5 have 
been emphasized. The colours of galactic stars have been super- 
imposed: light grey stars are those from the model data of Bessel, 
Castelli & Plez (1998), while the darker grey stars are from the 
photometric standards of Landolt (1992). Black stars represent 
M dwarfs and use data obtained from Bessell (1991) and Leggett 
(1992). 

Our selection criteria are indicated by the dashed lines. A broad 
area of colour-colour space is enclosed, allowing for efficient selec- 
tion of candidate high-redshift galaxies in spite of photometric er- 
rors, as long as a sufficiently large break in f7 — B can be observed. 
In practice this as the limiting factor, owing to the relatively poor 
quantum efficiency U of the TEK CCDs used compared to that 
of more modern detectors. 

Galaxies at redshifts above 2 3.5 will have almost complete 
absorption in U and thus extreme U — B colours off the top of 
this plot, but in practice we will only measure lower limits to this 
colour. As redshift increases, increasing absorption in B will push 
B — R further towards the red. Making the upper bound to the 
cut in B — R (or more generally for the two longer wavelength 
bandpasses) further to the red allows for selection of objects at 
higher redshifts than are shown on this plot, however this must 
not be pushed so far as to reach the stellar locus. 

6 MODEL COLOUR-COLOUR PLOTS 

Colour-colour plots for some of the three-filter photometric 
systems used in our observations are presented in Figures ^ 
to ^ A wider range of plots is available via the WWW at 
our website (Stevens 1999a), together with an interface for 
the generation of plots for arbitrary colour pairs within a 
substantial database of optical, NIR and HST filters. 

Galaxy colours are plotted at redshift increments of 
0.05, together with colours typical of galactic stars, and our 
chosen selection criteria are shown. Young galaxies within 
the targetted redshift ranges have been highlighted for clar- 
ity. 

In producing these plots, no assumptions have been 
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Figure 6. 

As Figure ^ but for U — G and G — R. G magnitudes for stars are 
derived by interpolation from the B and V magnitudes. Note in 
particular tlie decrease in the area of colour-colour space in the 
selection area relative to Figure B 
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Figure 7. 

As Figure |5| but ior U - V and V - R. The low V — R colours 
for stars and low-redshift galaxies mean that the colour must be 
tightly constrained to avoid severe problems with contamination 
of candidate high-redshift objects and use of these colours to 
select high-redshift galaxies therefore not recommended. However 
it may provide a useful check in deep multicolour surveys. 



Figure 8. 

As Figure ^ but for B ~ V and V — I, and with galaxies in the 
redshift range 3.5 < z < 4.5 emphasized. Again there is good 
separation between the high-redshift galaxy population and the 
low redshift galaxies and stars. 
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Figure 9. 

As Figure ^ but for B — R and R — I. In comparison with the 
B VI method (Figure ph , the separation between the high-redshift 
galaxies and the stars and low redshift galaxies is poor, making 
efficient selection harder without good limits on the B — R colour. 
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Figure 10. 

As Figure ^ but for V — R and R — I. Selection of candidates 
at z ~ 5 is perfectly feasible, but spectroscopic confirmation of 
objects becomes increasingly difficult as the redshift increases. 

made as to whether objects are hkely to be observable. 
Galaxies are shown at all redshift intervals for which the 
ages of the oldest stars are less than the age of the Universe 
in any reasonable cosmology. In some cases evolved high- 
redshift objects will therefore be shown which are unphys- 
ically old or at least exceedingly unlikely to be observable 
within the foreseeable future. 

Such plots can be used to determine the most effective 
use of broad-band filters to detect Lyman-break objects in 
future surveys. For efficient selection one must try to sepa- 
rate the locus of Lyman-break objects as much as possible 
from that of stars and low-redshift galaxies, an important 
factor is thus to ensure a reasonable span in wavelength be- 
tween the filters used for selection. Figure |^ clearly shows 
the problem in using V and R as the two longer-wavelength 
bandpasses, since the two are too close in effective wave- 
length to give a significant colour term on intrinsically red 
objects. Using G as the central filter is better (Figure H) but 
Figure ^ shows that B is optimum. 

Obviously the choice of filters will be influenced some- 
what by the targetted redshift, three-fllter selection working 
most effectively where the central filter spans the wavelength 
of redshifted Lya. UBR will be ineffective at selection of ob- 
jects much beyond z ~ 3.1, whereas Steidel and co-workers 
have already shown the effectiveness of a system similar to 
our UGR out to around z « 3.5. Again, UGR is inappropri- 
ate beyond this: for instance we note that the radio galaxy 
4C 41.17 aX z = 3.80 is too red in G — J? to be selected in 
UGR, but is a clear dropout in UVR. 

When searching for objects beyond z 3.5, it becomes 
preferable to use B rather than U as the shortest wavelength 
bandpass. While UVI selection would in principle work well, 
the increased observing time required to make sufficiently 



deep observations in / rather than R owing to the increased 
sky background makes it advantageous to use B in place 
of U to decrease the necessary observing time in the bluest 
filter. As a central wavelength, V is preferable to R owing 
to the greater separation between high-redshift galaxies and 
low redshift objects, as can be seen from Figures H and 0. 



7 EXTENSIONS TO Z > 4.5 

Selection of candidates at redshifts beyond z ~ 4.5 is per- 
fectly feasible given sufficiently deep observations, although 
spectroscopic confirmation of such objects becomes much 
more difficult due to absorption and emission lines becom- 
ing redshifted to among the OH bands of the near-IR. 

Steidel's group have added an I filter to their UnGTZ 
system in order to look for G dropouts at z ~ 4.5, and have 
had considerable success (Steidel et al. 1999). At higher red- 
shifts still, V dropouts can be searched for. A VRZ filter sys- 
tem is in theory very effective at selecting galaxies at z ~ 5, 
but observational considerations will in general make this 
impractical owing to the very high sky background, poor 
sensitivity of CCDs and severe problems with fringing if 
thinned CCDs are used. 

In practice VRI is likely to be more suitable, indeed we 
use the data obtained in one of our fields to conduct a pilot 
study, but conclude that our data are insufficiently deep 
to probe a population of V dropouts at z ~ 5. Z imaging 
could be more effective with the use of large-format infra-red 
arrays, which are a factor of ~ 2 more sensitive than CCDs 
at such wavelengths, but the much greater area of CCDs, 
particularly in the form of large mosaic cameras, will give 
them the edge for some time to come. 

Deep multicolour imaging surveys using large-format 
detectors will reveal many candidates. One such survey is 
being undertaken by researchers at Oxford using the Wide 
Field Camera on the Isaac Newton Telescope, imaging ten 
square degrees in BVRI to depths comparable with this 
work, with subareas also imaged in U and Z. Among the 
many possibilities of these data is that of searching for bright 
Lyman-break galaxies spanning the range 2.8 ^ z ^ 5. Given 
the detection limits of the survey, we expect to find ~ 15000 
U dropouts at z ~ 3, 2000 or more R dropouts at z ~ 4 and 
~ 100 V dropouts at 2 ~ 5 (Stevens 1999b). 

Shallower surveys can still make use of similar tech- 
niques for finding high-redshift objects. Though insuffi- 
ciently deep to find galaxies at such redshifts, the Sloan 
Digital Sky Survey will probe quasars and other bright ac- 
tive galactic nuclei. Very similar colour selection techniques 
to ours have been applied to commissioning data to iden- 
tify candidate quasars at redshifts above z }t 3.6 (Fan et al. 
1999; Fan et al. 2000). 

Beyond z ^ 6, intergalactic attenuation will affect even 
Z band and objects at z > 6 with star-formation rates typ- 
ical of objects at 2 ~ 3 will have NIR magnitudes of 26 or 
beyond, difficult to observe even with long observations on 
8-lOm class telescopes. Ultimately observations of the ear- 
liest star-formation in the Universe at 10 J; z ;S 30 must be 
the preserve of NGST and similar instruments. 
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x" IfcJlQ 


Z 


Filter choice 

Ideal of UBVRIZ 


? 


> 5.3 


r'i'z' 


RIZ 


TN J0924^2201 


5.19 


VRI 


VRI 


6C 0140+326 


4.41 


GRI 


BVI 


8C 1435+635 


4.25 


GRI 


BVI 


4C 41.17 


3.80 


BVI 


BVI 


6C 0032+112 


3.65 


BVI 


BVI 


B2 0902+343 


3.40 


UGR 


UBR 




^ 3.0 


u'GR 


UBR 



Table 2. Choice of filters for searching for galaxies at sim- 
ilaj: redshift to the central radio galaxy in various high— 
redshift radio galaxy fields, including those targetted in 
our own work (Stevens, Lacy & Rawlings 2000). u' g'r'i' z' 
refer to the SDSS filter set. At 2 ^ 3 a filter slightly bluer 
than the standard U filter such as u' or the Un of Steidel 
et al. will give a sharper break below Lya, while the lesser 
overlap between the filters of the Sloan set are more suit- 
able when searching for R dropouts. Inevitably the choice 
depends on the availability of filters for a particular in- 
strument, but in many cases use of a customized filter set, 
with minimal overlap between filter bandpasses, will pro- 
vide greater efficiency at a particular target redshift. In 
particular we conclude that such techniques are likely to 
perform very poorly for target redshifts 4 ^; z J; 4.7 if one 
is restricted to a standard set without a G or similar filter 
between B and V. 
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